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ABSTRACT 

, At a density near a few xlO^ g cm"^, the subsonic burning in a Type la 
supernova enters the distributed regime (high Karlovitz number). In this regime, 

1 ^ . turbulence disrupts the internal structure of the flame, and so the idea of 1am- 
pq . inar burning propagated by conduction is no longer valid. The nature of the 

burning in this distributed regime depends on the turbulent Damkohler number 
(DaT), which steadily declines from much greater than one to less that one as 
the density decreases to a few xlO^ g cm~^. Classical scaling arguments pre- 
. diet that the turbulent flame speed st, normalized by the turbulent intensity ii, 

Q_i. follows st/u = Da^. for Day -1- The flame in this regime is a single turbulently- 

Q ! broadened structure that moves at a steady speed, and has a width larger than 

+i ! the integral scale of the turbulence. The scaling is predicted to break down at 

ci '. Da-r ^ 1, and the flame burns as a turbulently-broadened effective unity Lewis 

number flame. This flame burns locally with speed sx and width /a, and we refer 
to this kind of flame as a A-flame. The burning becomes a collection of A-flames 
spread over a region approximately the size of the integral scale. While the total 
T^ ! burning rate continues to have a well-defined average, st ~ u, the burning is 

^ ! unsteady. We present a theoretical framework, supported by both ID and 3D 

^ [ numerical simulations, for the burning in these two regimes. Our results indicate 

that the average value of st can actually be roughly twice u for Dar^l, and 
that localized excursions to as much as five times u can occur. We also explore 
the properties of the individual flames, which could be sites for a transition to 
rS I detonation when Dar ~ 1- The A-flame speed and width can be predicted based 

5^ I on the turbulence in the star (specifically the energy dissipation rate e*) and the 

turbulent nuclear burning time scale of the fuel t^^^. We propose a practical 
method for measuring sx and Ix based on the scaling relations and small-scale 
computationally-inexpensive simulations. This suggests that a simple turbulent 
flame model can be easily constructed suitable for large-scale distributed super- 
novae flames. These results will be useful both for characterizing the deflagration 
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speed in larger full-star simulations, where the flame cannot be resolved, and for 
predicting when detonation occurs. 

Subject headings: supernovae: general — white dwarfs — hydrodynamics — 
nuclear reactions, nucleosynthesis, abundances — conduction — methods: nu- 
merical — turbulence — distributed flames 



INTRODUCTION 




In I Asp den et al.l ( l2008al ) (henceforth Paper I) three-dimensional simulations of resolved 
flames were performed that examined the interactions between turbulence and a carbon- 
burning flame in a Type la supernova at different densities. Because of the strong dependence 
of flame width and speed on density, this study can be viewed as a survey of the behavior 
of the flame at variable Karlovitz number, 

Ka = 

Here sl and l^ are the laminar flame speed and width, respectively, and u and / are the turbu- 
lent intensity (rms velocity) and integral length scale (deflned conventionally as the integral 
of the longitudinal correlation function). For Kasl, the flame is laminar with propagation 
determined by a balance between conductive and burning time scales. Such laminar flames 
have a large Lewis number, which is to say thermal diffusion occurs much faster than carbon 
diffusion. This means that perturbing the flame surface into the ash leads to a focusing of 
heat by diffusion, enhancing the burning rate, burning away the perturbation. Similarly, per- 
turbing the flame surface into the fuel leads to a defocusing of heat by diffusion, decreasing 
the burning rate. As a result, the flames are thermodiffusively stable. At small-to-moderate 
Karlovitz numbers (Ka^l), it was found that this thermodiffusively-stable nature led to a 
balance between local enhancement of the flame and local extinction, and so the turbulent 
flame speed remained close to the laminar value. 

However, once the Karlovitz number was sufficiently large (Ka^lO), turbulence was 
sufficiently strong that the turbulent mixing dominated thermal diffusion, and the flame 
became completely stirred (it resembled a turbulent mixing zone) and its width was greatly 
broadened. The turbulent flame width was also much greater than the integral length scale. 
While the local burning rate was greatly reduced, the overall flame speed was a factor of flve 
or six times the laminar flame speed due to the enhanced volume of the burning. 

A key result from Paper I was that at high Karlovitz number (Ka ~ 230), turbulence 
dominates the mixing of fuel and heat while thermal diffusion plays only a minor role. 
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Figure [T] shows a joint probability density function (JPDF) of fuel and temperature from 
Paper I. The solid red line shows the distribution from the fiat laminar flame where thermal 
diffusion is the dominant mixing process, and the solid black line shows the distribution of 
fuel burning isobarically with no thermal or species diffusion. The agreement with the JPDF 
demonstrates that the turbulent flame is burning at an effective Lewis number close to unity 
- the mixing is chiefly due to turbulence. 

The simulations in Paper I were performed in small domains to ensure that the laminar 
flame was fully-resolved - the domain size was approximately twenty-five times the laminar 
flame width. The aim of this paper is to investigate turbulent flame speeds in larger domains, 
with the aim of predicting high Karlovitz flame speeds in a full star, and specifically testing 
predictions made in Paper I. In this paper, we will refine these predictions by giving an 
in-depth theoretical description of the distributed burning regime, which are then compared 
with one- and three-dimensional simulations. 



THEORETICAL DESCRIPTION OF THE DISTRIBUTED BURNING 

REGIME 



There is diversity in the literature when naming the burning regimes for large Karlovitz 
number. Here we use the term "distributed burning regime" (or "distributed reaction zone" ) 
to refer collectively to all burning with large Karlovitz number (Ka^tlO). We also follow 
Woosley et al.l ( l2009l ) in further subdividing this region based upon the turbulent Damkohler 
number. Day, to be defined below. S pecifically, Par < 1 will be referred to as the "well- 
stirred reactor" (see 
fiame regime" (see 



ar 

ill 



KersteinI (120011 ) or lWooslevet al.l (120091 )1 



PetersI ( 1986 ) or Wooslev et al. (120091) ) and Bar > 1 as the "stirred 



2.1. The Well-Stirred Reactor 



Damkohlerl (jl940l ) identified two b urning r egime s , whic h he referred to as "small-scale" 
and "large-scale" turbulence, see also IPeterd (Il999l . I2OOOI ). In the large-scale turbulence 
regime, laminar fiames are dragged around by large turbulent eddies, w hich is an appropriat e 
description of the burning in a Type la supernovae at high density (IWoosley et al.l 120091 ). 
The high Karlovitz number case from Paper I is in the "small-scale" turbulence regime, 
where the small-scale mixing is dominated by turbulence, rather than diffusion. 

In this regime, the turbulent eddies are strong enough to disrupt the internal structure 
of the fiame. The time scale of these eddies is faster than the time scale of the fiame, so 
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the flame is mixed before it can burn. Tlierefore, tlie burning takes place on an inductive 
time scale, or turbulent nuclear time scale, which is much slower than the laminar nuclear 
time scale. Damkohler predicted (by analogy with laminar flames) that the turbulent flame 
speed St and width It should depend upon the turbulent diffusion coefficient and the nuclear 
burning time scale. 



St 




and It = V^tT^uc^ (2) 



where t^^^ is the turbulent nuclear time scale (rj^^. = It/st), and Vt is the turbulent diffusion 
coefficient Vt = otul (not to be confused with Day), where a is an order one constant. For 
convenience, take a = 1, which can be thought of as absorbing the constant into the definition 
of the turbulent flame width, which is ambiguous. 

Keeping the Karlovitz number fixed (so the energy dissipation rate of the turbulence 
e* = u^/l is constant) and assuming that r'^^^ is constant (i.e. the limiting value has been 
achieved), then there is only one free parameter. This parameter can be written as a turbulent 
Damkohler number, defined as 

Day = ^^ = — — . 3 

^nuc u It 

It is straightforward to show that Da^ oc P'^, where the constant of proportionality is 
{s*^ t'^uc)"^ ■ Therefore, the free parameter can be thought of in terms of the integral length 
scale. Both st and It are proportional to Vji , so writing T>t = ul = e*^'^/^/^, it follows 
immediately that both st and It scale with /^/^ (or equivalently Da^). 



2.2. The Stirred-Flame Regime 

When the turbulence time scale tt becomes comparable to the turbulent nuclear time 
scale Tj^j., i.e. Da^ ~ 1, the turbulence on the integral length scale can no longer mix the 
flame before it burns. Therefore, the flame can not be broadened any further and a limiting 
behavior is reached. Specifically, for Dar ^1, the flame burns like a unity Lewis number flame 
(on the scale of the flame width) with local flame speed sx and width Ix- We refer to this 
kind of burning as a A-flame. 

Defining a turbulent Karlovitz number as 



'n3/ 



KaT=W^f, (4) 



it can be shown that Da^-Ka^- = a = 1. Therefore, when the scaling relations break down 
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at Dar ~ 1, it follows imniediately that the turbulent flame speed is equal to the turbulent 
intensity and the turbulent flame width is equal to the integral length scale. 

It is the limit Da^ ~ 1 that divides the distributed regime into the well-stirred reactor 
regime and the stirred-fiame regime. In particular, note that it is the turbulent Damkohler 
number that is the divide, and that Da^ = crDa^, where a = Tj^ic/t^c is the ratio of the 
nuclear time scales. Therefore, a A-flame can only exist in the stirred-fiame regime, i.e. 
Da-r > 1 and Ka^lO. 



Figure [2] shows a regime diagram, based on iPetersI (119991 . 120001 ). where we emphasize 
the divide in the distributed burning regime. The diamonds denote the simulations from 
Paper I, and the squares denote the simulations from the present paper (to be defined below). 
The circle denotes the intersection of the Ka = 230 line with Da^- = 1, which denotes the 
A-point, where the turbulent intensity and integral length scale are equal to the turbulent 
flame speed (sa) and width (/a), respectively. 

We emphasize that the A-flame speed and width are local measures, i.e. the flame burns 
at Sa on a scale of /a- These quantities will also vary due to turbulent intermittency. The 
A-flame will respond to the local turbulence, specifically e on the scale of Ix. Note that due to 
this response, the turbulent Gibson scale Iq = s\/{v? /I) is equal to the A-flame width, /^ = l\, 
and the Karlovitz number based on the A-flame is always one, i.e. Ka^ = {u'^lx) / {slJ,) = 1. 



The overall turbulent flame speed will be greater than sa due to enhanced flame surface 
area, and will resemble Damkohler's large-scale turbulence re gime, d u e to the presence of 



multiple A- flames across an integral length scale. Following iPetersI (119991 ). for example, 
the following simple expression, based on an enhanced flame surface area, can be used to 
illustrate the scaling behavior in this regime, 

^ = 1 + -. (5) 

Sa sx 

Specifically, in the limit of high Damkohler number, the turbulent burning speed tends to 
the turbulent intensity, i.e. st ^ u as Day — ?> oo. 

In the next section, we present calculations to investigate the turbulent flame speed as 
a function of Damkohler number for a high Karlovitz number flame. Specifically, we are 
looking for a relation of the form 

? = <^ (DaT) , (6) 

u 

for some dimensionless function if. Three burning regimes are to be expected. First, for 
Da^^l, the turbulent flame speed and width are predicted to scale with 

-^ = Da^- , and — = Da^- . (7) 

(Jj t 
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Second, for Day ~ 1, the flame is predicted to reach a hmiting behavior with speed sx and 
width Ix (on the scale of Ix)- Finally, for Da^ ^ 1, the flame is predicted to burn as a 
A-flame, where the overall flame speed is expected to be a few times the turbulent intensity, 
tending towards it as Da-p tends to infinity. 



2.3. The Turbulent Nuclear Time Scale 

Previous work has i mplicit l y assu me d that t he nu clear time scale rema ins unaffected by 



the tu rbulence, again see iPetersI (119991 ) or lPetersI fl2000l ). This was the case in lRopke fc Hillebrandt 



(120051 ) . where the relation 



/^ X 1/2 / ., N 1/2 

St / Vt \ ( ul ^ ' 



was used to derive a turbulent flame speed for a level-set method. However, as discussed in 
Paper I, due to the different distributions of carbon and temperature, the nuclear time scale 
is around an order of magnitude longer in the turbulent case. Allowing for different nuclear 
time scales, equation (|8]) becomes 



9^ / f,/ \ 1/2 / L xl/2 



Sl \slIlJ \r. 



-T 
' nuc 



(9) 



In Paper 1, the turbulent nuclear time scale was estimated using It/st, based on an estimate 
of the turbulent flame width. It was shown that this reduces the estimate for the turbulent 
flame speed in equation [H] by a factor of approximately cr^/^ ^ 0.3. A more refined approach 
for estimating rj^^, is given below. 

A corollary to the relation Da^Ka^ = 1, is that the turbulent nuclear time scale can 
be derived from a single measurement in the well-stirred reactor regime, specifically, the 
turbulent flame speed and properties of the turbulence alone. In particular, if we consider 
a reference case that is easily computed (e.g. case (e) from Paper I), where the turbulent 
intensity is uq and integral length scale Iq, then only the turbulent flame speed s^ needs 
to be measured. Assuming that the limiting nuclear time scale has been achieved, then 
the relation Da^Ka^ = 1 means that the turbulent flame width is l^ = UqIq/s^, and the 
turbulent nuclear time scale is r^^^, = l^/s^ = uqIq/s^ . For case (e) in Paper I, this gives 
h ^ 190cm, rj„, ^ 0.0098s, and a ^ 0.067 (which gives a^/^ ^ 0.26). 
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SIMULATION DESCRIPTION 



3.1. Three-Dimensional Simulations 



As in Paper I, we use a low Mach number hy d rodyn ainics code, adapted to the study 
of thermonuclear flames, as described in lBell et al.l (J2004f ). The advantage of this method is 
that sound waves are flltered out analytically, so the time step is set by the the bulk fluid 
velocity and not the sound speed. This is an enormous efficiency gain for low speed flames. 
The input physics used in the present simulations is largely unchange d, with the exceptio n 
of the addition of Coulomb screening, taken from the Kepler code (jWeaver et al.l Il978l ). 
to the ^^C(^^C,7)^'*Mg reaction rate. This yields a small enhancement to the flame speed , 
and is included for completeness. The conductivities are those reported in iTimmesI ( l2000l ). 
an d the equation of s t ate is the Helmholtz free-energy based general stellar EOS described 



m iTimmes &: Swestyl (J2000l ). We note that we do not utilize the Coulomb corrections to 
the electron gas in the general EOS, as these are expected to be minor at the conditions 
considered. 

The basic discretization combines a symmetric operator-split treatment of chemistry and 
transport with a density-weighted approximate projection method. The projection method 
incorporates the equation of state by imposing a constraint on the velocity divergence. The 
resulting integration of the advective terms proceeds on the time scale of the relatively slow 
advective transport. Faster diffusion and chemistry processes are treated time-implicitly. 
This integration scheme is embedded in a parallel adaptive mesh reflnement algorithm 
framework based on a hierarchical system of rectangular grid patches. The complete in- 
tegration algorithm is second-order accurate in space and time, and discretely conserves 
species r nass and entha l py. T he details of the adaptive i ncompressibl e flow solver can be 
found in lAlmgren et al.l ( 119981 ) , the reacting flow solver in iDay fc Belli (|2000[ ) , extension to 
generalized equation of state in iBell et al.l (|2004) . and an application to the Rayleigh- Taylor 
instability in type la SNe in IZingale et al.l (120051 ). 



The non-oscillatory flnite-volume scheme employed here permits the use of implicit 
large eddy simulation (iLES). This technique captures the inviscid cascade of kinetic energy 
through the inertial range, while the numerical error acts in a way that emulates the dissipa- 
tive physical effects on the dynamics at the grid scale, without the expense of resolving the 
entire dissipation sub r ange. An overview of the technique can be found in ICrinstein et al. 
(120071 ). I Asp den et al.l (l2008bl ) presented a detailed study of the technique using the present 
numerical scheme, including a characterization that allowed for an effective viscosity to be 
derived. Thermal diffusion plays a signiflcant role in the flame dynamics, and so is explic- 
itly included in the model, whereas species diffusion is signiflcantly smaller, and so is not 



explicitly included. 

The turbulent ve locity field was maintained using the forcing term used in Paper I and 



Asp den et al.l (l2008bl ). Specifically, a forcing term was included in the momentum equations 
consisting of a superposition of long wavelength Fourier modes with random amplitudes 
and phases. The forcing term is scaled by density so that the forcing is somewhat reduced 
in the ash. This approach provides a way to embed the flame in a zero-mean turbulent 
background, mimicking the much larger inertial range that these flames would experience 
in a type la supernova, without the need to resolve the large-scale convective motions that 
drive the turbulent e nergy cascade. The effects of resolution were examined in detail in 



Asp den et al.l (l2008bl ). where it was demonstrated that the effective Kolmogorov length 
scale is approximately 0.28Ax, and the integral length scale is approximately a tenth of 
the domain width. This is particularly relevant to the present study, because turbulence is 
the dominant mixing process. This means that the ILES approach can be used to capture 
the effects of turbulent mixing, which occurs on length scales much larger than the actual 
Kolmogorov length scale in the star. 

Figure [3] shows the simulation setup. The simulations were initialized with carbon fuel 
in the lower part of the domain and magnesium ash in the upper, resulting in a downward 
propagating flame. A high-aspect ratio domain was used to allow the flame sufficient space 
to propagate. Periodic boundary conditions were prescribed laterally, along with a free-slip 
base, and outflow at the upper boundary. 



3.2. One-Dimensional Simulations Using the Linear Eddy Model 



Simulations were also run using the Linear Eddy Model (LEM) of iKersteid ( 119911 ). 
This approach simulates the evolution of scalar properties in a one-dimensional domain, 
which can be interpreted as a line of sight through a three-dimensional turbulent flow. 
Diffusive transport and chemical reactions are coupled in a model that represents the effects 
of real three-dimensional eddies through a so-called "triplet map" . The advantage of LEM 
is that much finer resolution and larger length scales can be explored inexpensively. This 
is particularly important for the present problem, where the range of length scales is very 
large and the degree of turbulence very high. LEM has been successfully applied to a large 
range of phenomena, especially turbulent terrestrial combustion. Its disadvantage is that 
when applied to a novel environment like a supernova, there are two overall normalization 
factors that must be adjusted, the effective turbulent dissipation rate and the integral scale. 
By using LEM for the present problem, we hope to satisfy two goals - first to show that 
very similar answers are obtained using quite different techniques and second, to calibrate 
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the uncertain constants in LEM for the supernova problem. 



LEM has previously been compared with the same three-dimensional code in lWoosley et al 



(120091 ). When the nuclear physics and fuel temperature were adjusted to be the same, it 
was found that best agreement occurred for an LEM constant C = 11 and an integral scale 
in LEM that was three times that in the numerical simulation. Those same values are used 
here except that C has been assumed to be 10. As expected, this prescription gives excellent 
agreement in the well-stirred reactor regime where it was calibrated, but, as we shall see 
underestimates the flame speed by as much as a factor of two in the stirred-fiame regime. 
Thus a normalization that depends on Day is appropriate, and C = 3 to 5 for Day > 1. 



4. RESULTS 

4.1. Three-Dimensional Results 

The high Karlovitz case from Paper I was used as the starting point for the present study, 
and is referred to here as case (a). Because turbulent diffusion dominates the mixing of fuel 
and temperature, the resolution requirements are significantly relaxed - thermal diffusion has 
to be resolved to capture the laminar flame correctly, but the turbulent diffusion coefficient is 
much larger, and so fewer cells are required to resolve it. In fact, keeping e* constant means 
that the diffusion coefficient scales with T>t ~ l^'^, and a mixing length can be defined 
analogous to the Kolmogorov length scale rjx) = {'D^/e*y^'^, which scales as rjx) ~ /• Thus, 
the smallest scales that need to be resolved become larger as the domain size grows. 

The approach taken to achieve larger Damkohler numbers (i.e. larger length scales) was 
to start with case (a), which has a domain width of 1.5 x lO^cm, and run the same case at 
an eighth of the resolution, i.e. with a low-resolution cell width of Axlr = Axhr/8. The 
turbulent flame speed was used as the primary diagnostic. Agreement in the turbulent flame 
speed provides confldence in capturing the effects of turbulent mixing (relying on the ILES 
approach) with the new cell width; it is the turbulent flame speed that is of primary interest. 
A higher Damkohler number was then achieved by running in a domain eight times larger 
with the new cell width, i.e. with a domain size of 1.2 x lO^cm. The turbulent intensity was 
adjusted accordingly to keep e* = u^/l constant (i.e. an eight-fold increase in length scale 
corresponds to a two- fold increase in velocity). The whole process was then repeated several 
times to reach a domain size of 6.14 x lO^cm, and span a range of Damkohler numbers from 
the base case of Da^ ?a 0.0064 to Da^- ^ 1.68, corresponding to cases (a) through (e). 

For larger length scales, the above argument no longer applies because the diffusion 
coefficient responsible for mixing fuel and temperature is approximately sxlx and does not 
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scale with ul for Da^^l; even case (e) is questionably resolved. However, we have included 
simulations (f) and (g) at Damkohler numbers 6.6 and 26, respectively, as an indication of 
what we can expect for DaT^^il. The simulation properties are summarized in tabled] 

Figure H] shows slices of density (top) and fuel consumption rate (bottom) for the seven 
cases (a)-(g) left-to-right. Note that the domain size increases by a factor of 8 each time. 
All of the figures have been normalized by the same values. Case (a) presents an extremely 
broad mixing region (much broader than the integral length scale), and the burning can 
be seen to occur at the high temperature (low density) end of the mixing zone. As the 
Damkohler number is increased, the relative width of the flame brush decreases as expected. 
Although the width of the flame brush appears to decrease, it actually increases, just more 
slowly than the domain size. For Dar^l, there appears to be a sharp interface between the 
fuel and products. This is because of the underresolved nature of the flames - the actual 
flame will be a broad mixing zone, but is not fully-captured here. 

Figure [5] shows the turbulent flame speed evolution for the seven cases. The flame speed 
is evaluated by finding the rate of change of the total fuel mass in the domain divided by 
the product of the cross-sectional area of the domain and the fuel density. 

The time scale has been normalized by the eddy turnover time tt = l/u. The solid lines 
denote the simulations at the full resolution (256^ cells in cross-section), and the dashed 
lines denote the simulations at the low resolution (32^ cells in cross-section). In most cases 
the low resolution simulations are in good agreement with the high resolution simulations 
(the only real outlier is case (d) where the low resolution simulation over-predicts the flame 
speed by 48%). Table |2] shows the mean flame speeds at both resolutions in each case. 

Figure |6] shows the turbulent burning speeds normalized by the turbulent intensity as a 
function of Damkohler number. Case (a) is the low Damkohler number on the left, increasing 
to the right to case (g). The marker denotes the mean in each case, and the vertical line 
denotes the range of values obtained during the averaging period (note that it is not the 
standard deviation). The solid line is the expected scaling relation from equation [71 The 
dashed lines denote Da-r = 1 and st = u, where the scaling relation is predicted to break 
down. Cases (a)-(e) are in very good agreement with the predicted scaling. We note that 
for cases (d) and (e), the flame speed measured is enhanced by the area of the flame that is 
burning, which can be increased by turbulence here because the flame brush is now smaller 
than the domain width; this explains the higher speeds obtained. The flames in these cases 
are too convoluted to extract a flame area to normalize by. Despite being significantly 
underresolved, cases (f) and (g) appear to show the expected break-down of the Damkohler 
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scaling; for Dar^l, the turbulence cannot enhance the flame speed according to equation 
([7]), and the normalized flame speed ceases to grow. 

Figure [7] compares the volumetric rate of burning for cases (f) and (g). Because these 
two cases are much less convoluted than cases (a)-(e) it was possible to extract an isosurface 
based on a temperature of lO^K. The solid line denotes the measured turbulent flame speed 
multiplied by the cross-sectional area of the domain, the dashed line denotes the measured 
flame surface area times the predicted burning speed sx. It appears that the flame speed is 
over-predicted in case (f) and under-predicted in case (g); the dash-dotted line denotes this 
speed multiplied by a factor of 0.67 for case (f) and 2.16 for case (g). We speculate that the 
disagreement is due to the underresolved nature of these cases, but maintain that s\ is still 
a reasonable estimate for a turbulent flame model. 



4.2. One-Dimensional Results 

LEM was used to simulate the same conditions as in Tabled] except that in each case the 
integral scale was multiplied by three and flner resolution was employed (Table [3]). Roughly 
40,000 time points were sampled for each case. The average flame speeds are plotted in Figure 
El For Dar ~1, the agreement with the three-dimensional simulations is excellent. Both show 
the same scaling relation as well as agreeing on the actual value of the speed. However, at 
large values of Da^n, the LEM results are almost a factor of two smaller. There could be 
several reasons for this. First, there are fundamental differences in the two approaches. 
LEM does not capture all of the multi-dimensional effects and has not been calibrated for 
this regime. Second, the resolution of the three-dimensional study is low and not all burning 
structures are well resolved. Studies with LEM suggest a mild dependence of st on the 
resolution with slower speeds at higher resolution. 



4.3. Variability of the Turbulent Flame Speed 

Figure Elcompares the probability density functions (PDFs) of the normalized local burn- 
ing rate in the three-dimensional and LEM calculations. To compare the three-dimensional 
results as closely as possible with LEM, the PDFs were evaluated by integrating the fuel 
consumption along line-outs, i.e. over individual columns of data as a function of height. 
Each PDF was then normalized by the mean burning speed. The PDFs have been shifted 
along the y-axis for clarity. In case (a), the PDF has a narrow Gaussian distribution cen- 
tered around the mean flame speed. As the Damkohler number increases, the PDF becomes 
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broader, i.e. a greater range of burning rates is observed locally. 

Both studies show a strong dependence of the spread of the PDF on Damkohler number. 
In the well-stirred reactor (Da^ < 1), the integral scale is less than the flame width. Many 
eddies turn over on the largest scale as burning moves through a region. The relation between 
temperature, carbon mass fraction, and location is smooth and the burning time well defined. 
There is only one flame structure. Consequently the speed does not vary greatly from the 
average. 

In the stirred flame, however, there are multiple regions of burning. Sometimes the 
burning is very fast, especially when a large single eddy envelops a new region of fuel. At 
other times it, almost goes out. The PDF thus shows a large spread in speed. Occasionally, 
the overall burning proceeds at a rate that is 2.5 to 3 times the average. Since the average 
itself if roughly twice u, this implies an overall burning rate up to six times u. 



5. CONCLUSIONS 

New one- and three-dimensional simulations have been presented to clarify and quantify 
the nature of carbon burning in a Type la supernova in the distributed burning regime. The 
characteristics of distributed burning depend upon the Damkohler number, and a range of 
Day from 0.006 to 26 has been explored. For Da^l, the expected scaling relations were 
demonstrated, 

— = Da^; , and — = Da^, . (11) 

u I 

For Da^^al, these relations break down and the turbulent flame reaches a maximum 
local flame speed sx and width Ix, measured on the scale of Ix- This A-flame interacts with 
turbulent eddies with length scales between Ix and the integral length scale I, which enhances 
the flame surfa ce area. The average turbulent flame properties are constrained by the large 



scales (see also lWoosley et al.ll2009[ ). This leads to an average overall turbulent flame speed, 
which to factor-of-two accuracy is given by the turbulent intensity, u. For the range of Da^l 
studied here, a good approximation is st = ^it. 



5.1. Consequences for Large-Scale Simulations 

Large-scale simulations (i.e. of a full-star) will require a turbulent flame model. One 
of the consequences of the results presented here is that the A-flame is ideally suited to the 
level-set approach, as it burns locally with speed sx- This means that a turbulent flame speed 
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can be evaluated in the following manner. Assume that the turbulence is sufficiently large 
and intense, i.e. in the stirred-flame regime (Da^ > 1 and Ka;jlO). Then the turbulence can 
be characterized by the energy dissipation rate e* = u^/l, and the burning depends on the 
turbulent nuclear burning time scale t^^^. The nuclear burning time scale is equal to Ix/sx 
and the turbulent intensity at the length scale of the A-flame is sx and so e* = s\/lx- Solving 
for Sx and Ix gives 



5a 



Ve*^, and h = V£*rJ„c'^- (12) 



Note these expressions do not require a fixed Karlovitz number (it's incorporated by e*), and 
the the length scale is the same as predicted in Paper I. To implement a level-set approach for 
a A-flame, Tj^^, should be evaluate d as a function of fue l density and temperature (using the 



approach outlined in section |2TT] or I Woosley et al.l ( l2009l )) and coupled with the (local) energy 



dissipation rate to evaluate the A-fiame speed according to equation ( !T2|) . If the resolution 
requirements of Paper I were applied to the A-fiame recovered in the present paper (i.e. 
Ix ~ 4Ax with 256^ cells in cross-section), an integral length scale could be achieved that is 
of order thousands of times larger than the original study. 



5.2. Consequences for Detonation 

The instantaneous flame speed for Da^ ^1 is irregular, with frequent excursions up to 2.5 
times the average, which is approximately twice the turbulent intensity for the simulations 
presented here. These variations could be crucial if a transition to detonation is to occur. 
Detonation will not happen in the laminar regime, i.e. at high density, because each fiamelet 
burns at a steady rate and has a thickness that is far below the critical mass for detonation. 
It is also unlikely to happen in the well-stirred reactor regime, because the burning is steady 
and slower than the turbulence on the integral scale, which in turn is subsonic. The burning 
time scale is very long in the well-stirred regime, making it very difficult to burn, for example, 
a 10 km region on a sound crossing time. 

Thus if detonation is to occur spontaneously, it is likely to happen in the stirred-flame 
regime. Even there though. Day should not be too large, or the mixed regions will not burn 
supersonically (recall st ^ u for large Da-r). The best conditions therefore occ ur where 



there are multiple A-flame s across an integral length scale, i.e. for small Da^ ^1 (jWoosley 



20071 : IWoosley et al.l |2009[ ) . The calibration of LEM determined by comparison with the 



three-dimensional studies suggest a normalization c onsta nt of C = 5 is most appropriate for 



Da > 1 and this is the value used by lWoosley et al.l (120091 ). This strengthens their conclusion 



that a detonation is possible if u exceeds about one-fifth sonic. 
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Fig. 1. — Joint probability density function of fuel and temperature from case (e) in Paper I. 
The solid red line shows the distribution from the flat laminar flame where thermal diffusion 
is the dominant mixing process. The solid black line shows the distribution of fuel burning 
isobarically with no thermal or species diffusion. 
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Fig. 2. — Regime diagram based on iPeterd ( 119991 . 120001 ) . Here we emphasize the separation 
of the distributed burning regime into the well-stirred reactor and stirred-flame regimes by 
the turbulence Damkohler number Day = 1- (LF - Laminar Flames, WF - Wrinkled Flames, 
CF - Corrugated Flames, TRZ - Thin Reaction Zone, WSR - Well-Stirred Reactor, SF - 
Stirred Flame). Together the WSR and SF regimes make up the distributed burning regime. 
The diamonds denote the simulations from Paper I, and the squares denote the simulations 
from the present paper. The circle denotes the intersect of the Ka = 230 line with Da^- = 1, 
which denotes the A-point, where the turbulent intensity and integral length scale are equal 
to the turbulent flame speed (sa) and width (/a), respectively. 
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Fig. 3. — Diagram of the simulation setup (shown in two- dimensions for clarity). The domain 
is initialized with a turbulent flow and a flame is introduced into the domain, oriented to that 
the flame propagates toward the lower boundary. The turbulence is maintained by adding a 
forcing term to the momentum equations. The top and bottom boundaries are outflow and 
solid wall boundaries, respectively. The side boundaries are periodic. 
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Fig. 4. — Two-dimensional slices showing density (top) and fuel consumption rate (bottom) 
for all cases (a)-(g) (left-to-right). All of the figures have been normalized by the same values. 
It is important to note the domain increases by a factor of 8 each time. In particular, the 
domain size in the final case is over two hundred and fifty thousand times larger than the 
first case; despite the burning rate looking reduced, it is actually greatly increased. 
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Fig. 5. — Turbulent flame speeds st as a function of time (normalized by turbulent time scale 
Tx- The solid lines denote the simulations at the full resolution (256^ cells in cross-section), 
and the dashed lines denote the simulations at the low resolution (32^ cells in cross-section). 
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Fig. 6. — Turbulent flame speeds normalized by turbulent intensity st/u as a function 
of turbulent Damkohler number Da-r. The solid black line denotes the expected scaling 
behavior for Da.T^l, bounded by the dashed lines. The flame speeds appear to be in good 
agreement with the predicted scaling, and despite the lack of resolution, appear to roll over 
when DaT^l. The crosses denote the mean speeds from the LEM calculations. 
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Fig. 7. — Volumetric rate of fuel consumption. The solid line is the measured rate Ast, 
where A is the cross- sectional area, the dotted line is the product of the limiting flame speed 
sx with the measured flame surface area, and the dash-dotted line is the latter adjusted by 
a factor of 0.67 for case (f) and 2.16 for case (g). There is a definite correlation between 
the flame surface area and the fuel consumption, but the estimate sx appears to be an over- 
prediction for case (f) and an under-prediction for case (g). We speculate that this is due 
to the underresolved nature of these cases and that sx is still a reasonable estimate for a 
turbulent flame model. 
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Fig. 8. — Probability density functions of normalized local burning rate for both the 3D 
(solid) and LEM (dashed) calculations. Each case has been shifted along the y-axis for 
clarity. 



Case 


(a) 


(b) 


(c) 


(d) 


(e) 


(f) 


(g) 


Domain width, L (cm) 


1.5 X 10^ 


1.2 X 10^ 


9.6 X 10^ 


7.68 X 10^ 


6.14 X 10^ 


4.92 X 10'^ 


3.93 X 10^ 


Domain height, H (cm) 


1.2 X 10^ 


4.8 X 10^ 


3.84 X 10^ 


3.07 X 10^ 


2.46 X 10^ 


1.97 X 10^ 


1.57 X 10^ 


Integral length scale, / (cm) 


1.5 X 10^ 


1.2 X 10^ 


9.6 X 10^ 


7.68 X 10^ 


6.14 X 10^ 


4.92 X 10^ 


3.93 X 10^ 


Turbulent intensity, u (cm/s) 


2.47 X 10^ 


4.93 X 10^ 


9.86 X 10^ 


1.97 X 10^ 


3.95 X 10^ 


7.89 X lO'' 


1.58 X 10^ 


Laminar Damkohler, Da^ 


9.32 X 10-2 


3.73 X 10-1 


1.49 X 10° 


5.97 X 10° 


2.39 X lOi 


9.54 X lOi 


^.82 X 102 


Turbulent Damkohler, Day 


6.23 X 10-3 


2.49 X 10-2 


9.97 X 10-1 


3.99 X 10-1 


1.60 X 10° 


6.38 X 10° 


2.55 X lOi 


High resolution (A^ = 256) Axhr (cm) 


5.86 X 10-1 


4.69 X 10° 


3.75 X lOi 


3.00 X 102 


2.40 X 10^ 


1.92 X 10^ 


1.54 X 10^ 


Low resolution (A^ = 32) Axlr (cm) 


4.69 X 10° 


3.75 X 10^ 


3.00 X 102 


2.40 X 10=^ 


1.92 X 10^ 


1.54 X 10^ 


1.23 X 10^ 



Table 1: Simulation properties. 



Case 


(a) 


(b) 


(c) 


(d) 


(e) 


(f) 


(g) 


High resolution (cni/s) 


1.84 X 10^ 


8.10 X 10^ 


3.51 X 10^ 


1.56 X 10^ 


5.94 X lO*' 


1.66 X 10^ 


3.49 X 10^ 


Low resolution (cm/s) 


1.86 X 10^ 


8.57 X 10^ 


3.91 X 10^ 


2.31 X 10^ 


6.5 X 10^ 


1.62 X 10^ 


3.16 X 10^ 


Percentage error 


1.1 


5.8 


11.4 


48.1 


9.4 


-2.4 


-9.5 



K2 



Table 2: Measured turbulent flame speeds. 



Case 


(a) 


(b) 


(c) 


(d) 


(e) 


(f) 


(g) 


Turbulent intensity, u (cni/s) 


2.47 X 10^ 


4.93 X 10^ 


9.86 X 10^ 


1.97 X 10^ 


3.95 X 10^ 


7.89 X 10^ 


1.58 X 10^ 


Integral length scale, / (cm) 


4.50 X 10^ 


3.60 X 102 


2.88 X 10^ 


2.30 X 10^ 


1.84 X 10^ 


1.48 X 10^ 


1.18 X 10^ 


Resolution, Ax (cm) 


4.89 X 10^1 


2.93 X 10° 


1.47 X 10^ 


9.77 X 10^ 


7.32 X 10^ 


2.44 X 10^ 


2.44 X 10^ 


Average speed (cm/s) 


1.96 X 10^ 


7.26 X 10^ 


2.50 X 10^ 


8.70 X 10^ 


2.85 X 10^ 


8.36 X 10^ 


2.14 X 10^ 



to 



Table 3: Characteristics of LEM studies. 



